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ABSTRACT 

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares 
to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been 
reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor 
imaged because of the instrumental limitations of previous radio telescopes. 

Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of 
S bursts, as well as their origin and possible emission mechanisms. 

Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales 
and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is 
currently limited to one image per second. 

Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups 
of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz s -1 and a wide range of circular 
polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic 
field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system 
that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. 

Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed 
at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the 
characteristics of electron-cyclotron maser emission. 
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1. Introduction 

The Sun is an active star that produces numerous large scale 
energetic phenomena often accompanied by solar radio bursts 
such as the Type I, II, III, IV and V radio bursts. At decime¬ 
tre/metre wavelengths, numerous other radio bursts have been 
observed, for e xample metric spikes (|Benz et al.||1982| |1996| ), 
drifting spikes (Elgaroy & Sveen 1919) and supershort solar ra¬ 
dio bursts ( Magdale nic et al.|2006 ). These bursts exhibit a variety 
of fine frequency structures and time scales of <1 s, which can 
be an indicator of small scale processes occurring in the solar 
corona. However, at longer wavelengths, corresponding to fre¬ 
quencies <100 MHz, there have been very few studies of short 
fine structure radio bursts due to sensitivity and imaging limita¬ 
tions of previous radio telescopes. 

At low frequencies, below 100 MHz, Ellis ( 1969| ) used 32 
broad-band dipoles to identify a new type of short duration ra¬ 
dio burst that he named fast drift storm bursts. Later, [McConnell) 
(1982]) observed these bursts using the Llanheme radio telescope 
consisting of 4096 dipoles ( [Ferris et aL||1980| ). He showed that 
their drift rate is slower than that of Type III radio bursts (about 


1/3 the drift rate of Type Ills). |McConnell| ( l982| ) renamed fast 
drift storms as “solar S bursts" due to their similarity to jovian S 
bursts. S bursts appear as narrow drifting lines in dynamic spec¬ 
tra with a temporal width of up to a few tens of milliseconds 
and total duration of about 1 s. They h ave been observed in a 
frequency range of 30-150 MHz ( |Ellis| 1982] |McConnell| 1982| ). 
They have short instantaneous bandwidths of about 120 kHz and 


instantaneous duration of 50 ms at frequencies of 40 MHz (Mc- 
Connell|1983| ). More recently, [Briand et ak] ( |2008| ) and |Melnik 


et al.| ( [2010[ ) have reported observations of S bursts in the fre¬ 
quency range 10-30 MHz. These were observed with the UTR-2 
radio telescope co nsisting of over 2000 dipoles with a tot al area 
of over 30000 m 2 ( Abranin et al. 2001 Melnik et al. 201 1) . Most 
S bursts have negative drift rates, i.e. from high to l ow freq uen- 
cies, however a few S bursts have positive drift rates (M cConnell] 
1982] |Briand et al.|2008| ). All S bursts occurred during times of 
other solar activity such as Type III and Type Mb radio bursts 


( |McConneII||1982[ |Briand et~aL]|2008[ |Melnik et al.||2010| ). The 

majority of S burst studies consider plasma emission as the most 
probable emission mechanism ( [Zaitsev & Zlotnik[1986t Melnik 


|et al.|[20T0| ). To date, S bursts have only been observed with 
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Fig. 1. LOFAR tied-array beam observations of Type III radio bursts and solar S bursts, (a) Map of 170 tied-array beams covering a field-of-view 
of ~ 1.3° about the Sun. The full-width half-maximum (FWHM) of the beams at zenith at a frequency of 60 MHz is represented by the red circle 
and the size of the optical Sun is represented by the yellow circle, (b) Dynamic spectrum recorded for a period of 1.3 minutes corresponding to 
Beam 14 containing Type III radio bursts, (c) Dynamic spectrum recorded for a period of 1.3 minutes corresponding to Beam 35 containing Type 
III radio bursts and S bursts. The two arrows indicate the beams that recorded the two dynamic spectra in (b) and (c) pointing at the time and 
frequency corresponding to the intensity values in (a). 
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Fig. 2. S -burst property distributions: start and end frequency (a), total 
frequency bandwidth (b), duration (c) and drift rate (d). The most fre¬ 
quently occurring value (mode) and the standard deviation are shown. 
The standard deviation was calculated using a normal distribution fit to 
each of the distributions. 


Table 1. S bursts characteristics estimated from LOFAR dynamic spec¬ 
tra are compared to results by |McConnell| ( |1982| . 



LOFAR 

McConnell 

Frequency range (MHz) 

20-80 

30-82 

Duration (s) 

0.5+0.3 

1.2 

Instantaneous duration (ms) 

20-400 

50 

Total bandwidth (MHz) 

2.5+1.3 

<5 

Instantaneous bandwidth (MHz) 

0.1-1.5 

0.1 

Drift rate (MHz s -1 ) 

-3.5+2.0 

-3 


to observe solar S bursts in high time and frequency resolution 
dynamic spectra and for the first time obtain spatial characteris¬ 
tics of S bursts from imaging. In this paper, LOFAR tied-array 
beams were used to study the spatial and spectral characteristics 
of S bursts. In Section 2 we give an overview of the LOFAR 
instrument and the observational method. In Section 3 we 
present the results of the spectral and imagining analysis which 
are discussed in Section 4. 


2. Observations and data analysis 


highly sensitive instruments as they have a low intensity com¬ 
pared to the majority of radio bursts. So far, there has been no 
imaging of S bursts for more detailed studies of their origins and 
emission mechanisms. 


With the recent development of the Low Frequency Array 
(LOFAR, van Haarlem et al. 2013] ), we now have the possibility 


On 9 July 2013 starting at 07:00 UT, over 3000 S bursts were ob¬ 
served during a LOFAR observation campaign that lasted until 
14:30 UT. These observations were made at times of medium so¬ 
lar activity in the form of three C-class flares and the presence of 
four /3yd active regions (NOAA 11784,11785, 11787, 11789) on 
the visible side of the solar disc. In this time interval, a long-lived 
noise storm was observed above active region NOAA 11785 at 
frequencies greater than 150 MHz. 
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Fig. 3. Location of S bursts at two separate frequencies, (a) LOFAR 4-second dynamic spectrum showing negatively and positively drifting S 
bursts. The white line marks the timestamp in the images below, (b) SDO/AIA 171 A image at 11:25:11 UT and the overlaid LOFAR source of the 
negative drifting S burst at a frequency of 64-65 MHz at the time marked by the white line in (a), (c) SDO/AIA 171 A image 11:25:11 UT and the 
overlaid LOFAR source of the positive drifting S bursts at a frequency of 51-52 MHz at the time marked by the white line in (a). An animation of 
the temporal evolution is available in the online edition. 


LOFAR represents a new milestone in low radio frequency 
arrays and was constructed by the Netherlands Institute for 
Radio Astronomy (ASTRON). It consists of -7000 antennas: 
Low Band Antennas (LBAs) operating at frequencies of 10- 
90 MHz and High Band Antennas (HBAs) which operate at 
110-240 MHz ( [van Haarlem et al.1|2013| ). These antennas are 
distributed in 24 core stations and 14 remote stations across the 
Netherlands and 9 international stations across Europe. 


In this paper, we used one of LOFAR’s beam formed modes 
(Sta ppers et al.|[20lH |van Haarlem et al.] 2013| ) in the LB A 
frequency range. Figure la shows 170 simultaneous beams 
which were used to observe the Sun covering a field-of-view 
of 1.3° centred on the Sun. The 24 LOFAR core stations were 
used to produce these beams. Each beam produces a high time 
and frequency resolution dynamic spectra (~10 ms; 12.5 kHz) 
as can be seen in Figures lb and lc, in which the radio bursts 
can be identified. Since each beam has a spatial location, it can 
be used to extract the intensity of the radiation at that specific 
location. These intensity values can be plotted as ‘macro-pixels’ 
onto the tied array map for a chosen time and frequency as seen 


in Figure la. The arrows in Figure 1 denote the time, frequency 
and beam where the intensity value was extracted. 


Radio sources can be imaged by interpolating between the 
intensity of radiation from all beams and producing tied-array 
images of radio bursts (for more details, see Mor osan et al. 
2014| . All images of S bursts were averaged over a frequency 
bin of 1 MHz to reduce spectral noise and over a time period of 
50 ms. Side lobes are often visible in tied-array images due to the 
coherent combination of tied-array beams, however, most side 
lobes have been excluded manually as they occur at the same 
time as the source and always at the same location some distance 
away from the source. Due to the size of tied-array beams there 
is an uncertainty in determining the exact location of the radio 
source as the source could be located anywhere within the beam. 
This uncertainty is estimated to be 0.2 R Q at 80 MHz and 0.6 R Q 
at 30 MHz and it corresponds to the half power beam width. The 
main source of uncertainty in beam pointing is the ionosphere, 
which is most pronounced during periods of elevated solar ac¬ 
tivity. During our observations, there was no significant X-ray 
activity other than a few C-class flares. In addition, the centroid 
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Fig. 4. S bursts source position and their relation to the coronal mag¬ 
netic field. The AIA 171 A image at 11:24:12 UT is overlaid PFSS 
magnetic field extrapolations. The white contours represent the LOFAR 
S bursts source at 11:24:56 UT and a frequency of 74-75 MHz. One 
of the trans-equatorial loops connecting the bipolar page to the active 
region is labelled in the figure. 


position of our sources was found not to vary significantly over 
many hours. 


3. Results 

3.1. S bursts spectral characteristics 

An example of S bursts that occurred during a time period of 1.3 
minutes on 9 July 2013 at 07:00 UT is shown in Figure lc. The 
S bursts were recorded simultaneously with other solar activity 
mainly in the form of Type III and Type Illb radio bursts, some 
of which did not originate at the same location as the S bursts as 
can be seen in Figures lb and lc. 

Over 3000 S bursts have been recorded during a time inter¬ 
val of 8 hours on the same day in the LB A frequency range. 
Due to a variety of shapes (various lengths, time profiles, fre¬ 
quency profiles) and various intensity levels, similar to Type 
Ill/Type Illb bursts, it was necessary to identify these bursts 
manually. Distributions of their characteristics such as the start 
and end frequency, total bandwidth, duration and drift rate are 
shown in Figure 2. These distributions are asymmetrical and 
therefore the mode of the distributions is shown instead of the 
mean due to the fact that they are skewed. Most S bursts had 
short durations of <1 s and instantaneous durations at a fixed 
frequency of <400 ms, some as small as 20 ms. They were ob¬ 
served at frequencies between 20 to 80 MHz and the majority 
had a total bandwidth of -2.5+1.3 MHz and instantaneous band- 
widths at a fixed time varying between 0.1-1.5 MHz. Very few S 
bursts had longer bandwidths >20 MHz. Their drift rates varied 
with frequency from -1 MHz s -1 at frequencies of 20 MHz to 
-7 MHz s -1 at 75 MHz (Figure 2) but most bursts had a drift rate 
of -3.5+2.0 MHz s -1 . The standard deviation of all distributions 
is a significant fraction of the mode, therefore there is a wide 
spread in duration, bandwidth and drift rate of S bursts. 
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Tab le 1 shows that these val ues are comparable to the find¬ 
ings of |McConnell (1982 1983| ). A few S bursts were found to 
have positive drift rates and an example is shown in Figure 3a. 
Positive drifting S bu r sts have also been reported before by|Mc- 


|Connell| ( |l982H1983| ); |Briand et al.| ( |2008] ); |Melnik et al.| ( |26lO| ). 

Full Stokes data were recorded starting at 11:40:00 UT for 
a duration of 1 hour. During this time several hundred S bursts 
were observed as well as Type III radio bursts. The observed data 
is uncalibrated and we cannot estimate the degree of polarisation 
accurately, however, we can compare the S bursts polarisation to 
the background Type III radio bursts. Type Ills emitted at the har¬ 
monic of the local plasma frequency are believed to be weakly 
polarised (<10%, |Dulk & Suzuki) 1980| . We have found a wide 
range of different circular polarisation degrees for S bursts from 
2-8 times more polarised than the accompanying Type III bursts. 
This agrees with previous studies in which the S bursts studied 
were found to have a strong component of circular polarisation 
( |McConnell| 1982| ). 


3.2. S bursts spatial characteristics 

An example of a number of S bursts analysed in more detail is 
shown in Figure 3. Figure 3 is a snapshot from the movie accom¬ 
panying this paper. Figure 3a shows a 4 second LOFAR dynamic 
spectrum containing S bursts with both negative and positive 
drift rates. Figures 3b and 3c show the location of these S bursts 
overlaid on extreme ultraviolet images (EUV) at 171 A from the 
Atmospheric Imaging Assembly (AIA; Lemen et al. |2012| ) on¬ 
board NASA’s Solar Dynamic Observatory (SDO). S bursts are 
imaged at frequencies of 64-65 MHz to sample the negatively 
drifting S bursts and 51-52 MHz to sample the positively drift¬ 
ing S bursts. 

Tied-array imaging can produce images at very high cadence 
(50 ms) which is necessary for radio bursts with durations of 
<1 s. However the beams have a FWHM of 12.6' at a frequency 
of 50 MHz as they can only be produced by the LOFAR core. 
Its 2 km baseline is too small to produce arcsecond beam sizes. 
As a result the source sizes appear large. We cannot detect any 
movement in the source centroids due to frequency change. Both 
negatively and positively drifting S bursts appear to originate at 
the same location. All S bursts analysed in this observation cam¬ 
paign on 9 July 2013 appear to originate at the same location 
with insignificant deviations in the source centroids (<0.1 R Q ). It 
is possible that this results from the reduced spatial resolution of 
tied-array beams imaging. 

We investigated the origin of these S bursts by relating the 
source locations to the coronal magnetic field. Figure 4 shows 
that the radio source does not seem to be directly located on 
top of known active regions but it appears to be close to the top 
of trans-equatorial loops connecting a complex active region in 
the southern hemisphere and a large area of bipolar plage in the 
northern hemisphere. The top of these equatorial loops is located 
at a distance of about -1.8 R Q from the solar centre. This dis¬ 
tance is calculated from the 3D PFSS extrapolation of the loop 
and it represents the radial height of the loop top from the sur¬ 
face. The spatial coincidence of the S bursts and the top of these 
equatorial loops in Figure 4 suggests that S bursts are associated 
with these loops. 


4. Discussion and conclusion 

Over 3000 S bursts were observed by LOFAR at times of 
medium solar activity and simultaneously with other radio emis- 
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sion such as Type III and Type Illb bursts. Most S bursts had 
very short durations (~0.5 s) and short bandwidths (-2.5 MHz). 
Their drift rates varied with frequency with the majority drifting 
at -3.5 MHz s -1 at frequencies of 20-80 MHz (Table 1) which is 
about 1/3 of the Type III drift rate. Type Ills drift at -11.33 MHz 


in the frequency range of 40-70 MHz according to Mann 


|& Klassen] ( [2002] ). The majority of S bursts are found at fre 
quencies between 60-70 MHz (Figure 2a), but this result may 
be influenced by the 60 MHz peak in the sensitivity of LBAs 
antennas ( [van Haarlem et al.]|2013| ). S bursts were found to be 
significantly more polarised than Type Ills. Calibrated polarisa¬ 
tion observations of S bursts will allow for these properties to be 
further studied in a follow on paper. 

For the first time, we were able to image solar S bursts. Tied- 
array imaging shows that all S bursts originate in the same spatial 
location with very little deviation of the source centroid. There 
seem to be no features in the solar corona at the location of S 
bursts sources other than long trans-equatorial loops connect¬ 
ing an active region in the southern hemisphere and an extended 
area of bipolar plage in the northern hemisphere. We have inves¬ 
tigated two possible emission mechanisms responsible for solar 
S bursts: plasma emission and electron-cyclotron maser (ECM) 
emission which are discussed below. 

So far, plasma emission has been proposed by Zaits ev & 


|Zlotnik| ( |1986| ) and |Melnik et al.| ( |2010|) as the mechan ism re¬ 
sponsible for emission of S bursts. |Zaitsev & Zlotnik| (Q986) 
suggest that electrons with velocities 10-20 times above their 
thermal velocity excite plasma waves near the upper hybrid reso¬ 
nance frequency. These waves are then scattered by ions produc¬ 
ing electromagnetic waves at the plasma frequency (f p ). Melnik 


[et al.| ( |2010| ) propose a model in which S bursts are generated 
by the coalescence of magnetosonic waves and Langmuir waves 
and derive a minimum magnetic field of 2 G necessary for the 
generation of these bursts at the plasma frequency (f p ). 

In our observation, S bursts appear to be associated with 
equatorial loops stretching to high altitudes (-1.8 R Q ) where 
plasma emission is likely the dominant emission mechanism 
( Gary & Hurford||l989| . Knowing the height of these loops, a 
density model can be used to find the electron density at this dis¬ 
tance. Using the electron density we can estimate the frequency 
of plasma radiation at the top of the loops using the following 
equation: 


/ = nf p = nC jN e . 


( 1 ) 


The emission frequency, /, is given by the local plasma fre¬ 
quency, f p , multiplied by the harmonic number n and is directly 
proportional to the square root of the electron density, N e , in 
cm -3 , where C = 8980 Hz cm 3/2 is the constant of proportional¬ 
ity. 

We used the radial electron density models of |Newkirk| 
( 11961[ ), [Saito et al.| ( |1977| ) and|Mann et al.|(|1999| ) and the time- 
dependent density model of |Zucca et al.| ( 2014| ) to estimate elec¬ 
tron densities at specific distances. These models predict electron 
densities in the range of 1.2xl0 6 -l.lxl0 7 cm -3 at -1.8 R Q . This 
corresponds to a frequency range of -10-30 MHz for funda¬ 
mental plasma emission and -20-60 MHz for harmonic plasma 
emission, frequency ranges in which S bursts are found. In addi¬ 
tion, S bursts show a drift rate dependence with frequency sim¬ 
ilar to that of Type III radio bursts |McConnell| f 1980| . A simi¬ 
lar drift rate dependence with frequency was also found in our 
analysis which will be discussed in a follow-on paper. It is there¬ 
fore possible that S bursts are a form of plasma emission. The 
high degree of polarisation of some of these bursts requires that 


they are emitted at the fundamental plasma frequency ( Wentzel 
[T9S4) . 

However, S bursts have a significantly slower drift rate than 
Type III radio bursts (about 1/3 of the Type III drift rate), very 
narrow bandwidths and short lifetimes and they are significantly 
more polarised. This is indicative of ECM emission. Jovian S 
bursts, which are similar in appearance to solar S bursts, orig¬ 
inate due to cyclotron maser emission in the flux tubes con¬ 
necting Io or Io’s wake to Jupiter ( [Hess et al. [20071 ). F° r ECM 
to occur in the solar corona, the electron-cyclotron frequency 
(fs = eB12nm e , where B is the magnetic field) has to be greater 
than the plasma frequency ( [Melrose] 199 1| ). Therefore we require 
a high B field at the emission site in the corona (> 20 G), as well 
as low electron plasma density (N e ). While the plasma density 
is low at large heights in the corona where S bursts occur, the B 
field is small (-0.25 G) as estimated from the models of |Zucca| 
[et al.| ( |2014| ). From our analysis, the location of these S bursts 
does not meet the condition for ECM to occur. Presently we are 
not able to draw a conclusion on the emission mechanism of S 
bursts and we hope that future imaging observations of S bursts 
will constrain these results. 

The tied-array beams prove promising in detecting millisec¬ 
ond duration radio bursts with very high frequency and temporal 
resolution. In addition to this, tied-array beam observations can 
show the positions of these radio bursts sources. The accuracy in 
these positions is only limited by the beam size. Ionospheric ef¬ 
fects were not significant during this observation. No movement 
was detected in the S bursts source position and the Quiet Sun 
position was constant during the observation. The spatial reso¬ 
lution of tied-array beams can only be increased by increasing 
the baseline of the LOFAR beam-formed modes (i.e., extend be¬ 
yond the LOFAR core) which would not have sufficient beams 
to cover the Sun and is also currently not available as an obser¬ 
vational mode. However, at such low frequencies, scattering in 
the corona may limit the resolution of the beam sizes. Another 
way to increase the spatial resolution when observing S bursts 
is to decrease the cadence of LOFAR’s interferometric snapshot 
imaging to <0.1 s which is also currently not ideal due to the 
high data volume rates. This, however, can be compensated by a 
coarser frequency resolution and it is worth testing in the future. 
Tied-array beam analysis has been previously applied to Type III 
radio bursts ( Mor osan et al.|2014| ) and it shows a high potential 
for studying S bursts and other radio bursts within the LOFAR 
spectral range. 
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